arXiv:1504.03726vl [astro-ph.HE] 14 Apr 2015 


The INTEGRAL view of A0535+26 during the giant outburst 

of February 2011. 


N. Sartore^, E. Jourdain^’^, and J. P. Roques^’^ 

CNRS - IRAP, 9 Av. Colonel Roche. BP 44^46,F-31028 Toulouse cedex 4, France 
Universite de Toulouse, UPS-OMP, IRAP, Toulouse, France 

nicola.sartoreSirap.omp.eu 


ABSTRACT 

A0535+26 is a slowly rotating pulsar accreting from the wind of a massive Be star, and that 
exhibits two cyclotron absorption lines in its X-ray spectrum, at about 45 and 100 keV, respec¬ 
tively. Unlike similar sources, no significant variations of the energy of its cyclotron lines with 
flux were observed to date. The bright outburst of February 2011 thus offers a unique occasion to 
probe this peculiar behavior at flux levels not yet observed with present-day instruments. Here 
we report on the spectral and timing analysis of the data from the spectrometer SPI on-board 
INTEGRAL collected during the outburst. At the peak of the outburst the estimated luminosity 
is ^ 4.9 X 10^^ergs“^. The fundamental cyclotron feature is detected at all flux levels, and its 
centroid energy is positively-correlated with the flux of the source, confirming that A0535-I-26 is 
accreting at a sub-critical regime. The correlation seems to fall off at ^ lO^^ergs”^, suggesting 
the transition from a Coulomb-stopping regime to a gas-mediated shock regime. From the timing 
analysis we found that the pulsar was spinning up during most of the outburst, and that the 
spin-up rate correlates with the flux of the source, albeit the correlation is steeper than the one 
expected from the standard disk accretion theory. Finally, we show that the pulse profile of 
the source changes dramatically as the flux increases. At high luminosity the profile is highly 
asymmetric, implying an asymmetry in the geometry of the accretion flow. 

Subject headings: X-rays: binaries; pulsars: individual (A0535-I-26) 


1. INTRODUCTION 


High-mass X-ray binaries (HMXB) are systems 
composed by a massive, early type star which pro¬ 
vides the matter to an accreting compact object, 
typically a highly-magnetized neutron star, B ^ 
10^^ — 10^^ G. This has been inferred from the 
presence in their X-ray spectra of absorption-like 
features, which are interpreted as cyclotron reso¬ 
nant scattering features (CRSF), where Eqrsf 
11.6 (R/IO^^ G) keV, and B is the magnetic field of 
neutron star. Such strong magnetic fields channel 
the accretion flow towards the magnetic poles of 
the neutron star. The resulting anisotropy induces 
a modulation in the radiation emitted by the neu¬ 
tron star, which is then observed as an X-ray pul- 
sar (see e.g. Wilmj2014 ; Revnivtsev fc Mereghetti 


2014 for recent reviews). 

The accreted matter may come either from 
an accretion disk, via Roche lobe overflow, or 
directly from the wind of the massive compan¬ 
ion. Wind-fed HMXBs are further sub-divided 
into two groups, super-giant and Be systems, de¬ 
pending on the properties of the donor star. In 
particular, Be stars are fast-rotating OB stars 
in which the outflow takes the form of a dense 
equatorial disk, as suggested by observations of 
broad emiss ion lines at optical/in frared wave¬ 
lengths (e.g. Porter fc Riviniu^l2003ll . Be/HMXB 
exhibit long periods of quiescence, interrupted by 
bursts of activity w hich are classified in two groups 
( Stella et al. 1986ll . The so-called type I, or nor¬ 
mal, outbursts occur quasi-periodically at orbital 
phases close to the periastron passage, and reach 
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a luminosity of ~ lO^^ergs”^ at most. They are 
likely related to an increase in the accretion rate 
due to the passage of the pulsar through the equa¬ 
torial disk of the companion. On the other hand, 
type II (giant) outbursts are less frequent, with 
no apparent correlation with the orbital phase, 
and can reach a luminosity close to the Eddington 
limi10, TEdd 1-38 x 10^® (M/Mq) ergs“^, where 
M is the mass of the neutron star. The origin of 
type II out bursts is not clear and still matter for 
debate (see iNakaiima et al. 2014 , and references 
therein). 

Interestingly, these changes in luminosity can 
be used to test the accretion theory on strongly- 
magnetized pulsars at different accretion rates. 
In accreting pulsars a radiation-dominated shock 
forms in the accre tion column at a lumin osity 
Lx ^ 10®®ergs“^ ( Basko fc SunvaevIllQTOh . As 
the matter passes through the shock, it gives away 
most of its kinetic energy, which is transferred to 
the radiation field. Below the shock, the accreted 
matter slowly sinks towards the surface of the pul¬ 
sar where it can be further decelerated. The fi¬ 
nal deceleration depends on the luminosity. For a 
certain critical value Lcrit, which depends on the 
properties of the neutron star, like e.g. the mass 
and t he strength of mag netic field (see Equation 
32 of iBecker et al.ll2012ri . the pressure exerted by 
the upward radiation flux is such that the accre¬ 
tion flow arrives at the surface of the pulsar with 
zero speed. Above this critical value the dynam¬ 
ics of the accreted matter is dominated by radia¬ 
tion pressure, and both the radiative shock and the 
line forming region are expected to move upwards 
where the magnetic held is lower. This in turn im¬ 
plies an anti-correlation between the energy of the 
cyclotron line and the luminosity of the source, 
as reported e.g. for the Be/HM XB systems 4U 
0115 - 1-63 and V0332-I-53 (see e.g. Nakajimael.^ 


the dec eleration occurs de creases with luminos¬ 
ity (see iBecker et al. 2012 , for details) and, be¬ 


20061 : TsvEfankov et al.ll2006l : iMowlavi et aLl2006ll . 

On the other hand, in sources below the criti¬ 
cal luminosity, the radiation pressure cannot stop 
completely the in-falling matter, which is instead 
decelerated by Coulomb interactio ns towards the 
base of the accretion column (Burnard et al 


19911 1 . In this case, the height of the region where 


^At this luminosity, the force exerted by the emerging radi¬ 
ation on the in-falling matter through Compton scattering 
balances the gravitational force. 


cause of the increase of magnetic field at lower 
altitudes, a positive correlation between the en¬ 
ergy of the cyclotron line and the flux is then 
expected. Such behavior has been observed e.g. 


Vela X-1 (IStaubert et al. 

20071 Klochkov et al. 

2012: DeCesar et al.l 20131: 

Fiirst et al. 20141). 


1.1. The case of A0535-|-26 

A0535-I-26 (A0535 hereafter) is a transient 
Be/HMXB discovere d by the Ariel V satellit e dur¬ 
ing a giant outburst ( Rosenberg et al.|l975ll. The 
syste ms lies at a distance d ^ 2kDC (ISteele et al.l 
1199811 from the Earth, and it is composed by 
a slowly rotating pulsar, Pspin 103 seconds, 
orbiting the BOIIIe star HDE 245770 in an ec¬ 
centric orbit (e = 0.47 ± 0.02), with a period 


Porb = 111 ± 0.1 days ( Finger et al. 2006ll . Af¬ 


ter its discovery, A0535 exhibited several bursting 
episodes of both types I and II. The hard X-ray 
spectrum of the source during these outbursts 
showed two prominent absorption lines at a bout 
45 and 100 keV (e.g. Caballero et al. 2013 . and 
references therein). 

At variance with the examples reported above, 
the energy of the fundamental cyclotron line in 
A0535 was found to be remarkably stable at ^ 
45keV in almost all observations of the source. 
The only exception occurred during a short flar¬ 
ing episo de in which a higher va lue, ^ 52 keV, was 
reported ( Caballero et al.l 2008l l. This increase in 
the observed cyclotron energy was coupled with 
changes of the pulse profile, and was explained 
by the occurrence of instabilities in the magne¬ 
tosphere of the pulsar leading to a sudden in¬ 
crease in the accretion rate and changes in the 


geom etry of the accr etion column (iPostnov et al 


2008). Acco r ding to Klochkov et al.l ( 2011 ). and 


Muller et al. ( 2013al) . A0535 exhibits a positive 


correlation of the cyclotron energy with the lu¬ 
minosity on a pulse phase-dependent basis. While 
this result suggests that the source is in the sub- 
critical regime, a definitive proof of the correlation 
of the cyclotron lines with the luminosity in A0535 
is still lacking. 

Dedicated observations of A0535 during a gi¬ 
ant outburst are thus important in order to extend 
the range of luminosity covered by spectroscopic 
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studies of the source. In this regard, the two previ¬ 
ous giant outbursts (May/June 2005 and Decem¬ 
ber 2009) were not well suited for the task due to 
Sun constraints, i.e. the source was too close to 
the Sun in the sky to be safely observed. On the 
other hand, the giant outburst of February 2011 
was well covered by several instruments, and it of¬ 
fers a unique possibility to reveal changes, or lack 
thereof, in the spectral properties of A0535 and 
their relation with the observed luminosity, thus 
probing the accretion regime of the source. An 
analysis of Swift/^kT and INTEGRAL/lSGRl 
data below 80 keV taken during th e outb urst has 
been presented bv ICaballero et all ( 2012h . where 
they reported that the energy of the fundamen¬ 
tal cyclotron line was stable on all but one of the 
INTEGRAL/ISGKl observations. 

In this paper, we report on the observations of 
the Be/HMBX A0535 during the giant outburst 
of February 2011, performed with the SPI instru¬ 
ment on-board INTEGRAL. In Section [2l we give 
a description of SPI, as well as of the observa¬ 
tions and the data reduction procedure used in 
the paper. In Section [H a pulse phase-averaged 
spectroscopic analysis is presented, completed by 
a timing analysis and a qualitative study of the 
energy-resolved pulse profiles. Finally, in Section 
m we discuss our results and present the conclu¬ 
sions. 


2. OBSERVATIONS AND DATA RE¬ 
DUCTION 


INTEGRAL (|Winkler et al.ll2003h is a mission 
of the European Space Agency devoted to the 
observation of the high-energy sky (from ^ keV 
to ~ MeV). The satellite was launched from 
Baikonur in October 2002 and put in a highly ec¬ 
centric orbit. The spectrometer SPI (SPectrom- 
eter on Integral), operating between 20 keV and 
8 MeV, is one of the two main instruments on¬ 
board INTEGRAL. A detailed description of the 
instrume nt and its in-flight pe rfor mances can be 
found in Vedrenne et al.l ( 2003ll and Rogues et al 


(120031). It consists of 19 high purity Germanium 
detectors maintained at a temperature of ^80 
K, which ensures an excellent energy resolution 
all along the allowed energy range. The instru¬ 
ment uses a tungsten-made coded-mask offering 
a 2.6° angular resolution over a 30° field of view 


(f.o.v.). Particle and photon background reduc¬ 
tion is achieved thanks to a 5 cm Germanate Bis¬ 
muth anti-coincidence shield. 


INTEGRAL observations are carried-out on 
a revolution-by-revolution basis, each revolution 
(rev. in short) lasting ~3 days, of which ^ 2.5 
days give useful data. Each observation is fur¬ 
ther sub-divided into shorter, ~2-4 kiloseconds 
exposures called “science windows” (sew). The 
low number of pixel s enforces the adop tion of a 
dithering procedure ( Jensen et al.ll2003ll . in order 


to better estimate the contributions of the back¬ 
ground and of the different sources in the f.o.v. 
The satellite aiming point is thus shifted by ^ 2° 
around the observed target at each sew, following 
a predefined pattern, typically a 5° x 5° rectangle. 


The giant outburst of February 2011 of A0535 
has been monitored in ToD (Target of Opportu¬ 
nity) mode by INTEGRAL in 8 consecutive rev¬ 
olutions (rev. 1021 to 1028), which caught es¬ 
sentially the maximum and the decline of the 
outburst. The total exposure of each observa¬ 
tion varied between ^40 and ^100 kiloseconds. 
In addition, the rise of the outburst was ob¬ 
served serendipitously during rev. 1019, a ^ 
186 kilosecond-long calibration observation of the 
nearby Grab Nebula (Figure [J). A log of the IN¬ 
TEGRAL observations of A0535 used in this paper 
is reported in Table [T] 


The spectrum of each source in the f.o.v. of 
SPI is evaluat ed by means of a sky mod el fitting 
procedure ('see ljourdain fc Rogues 20091 for a de¬ 
tailed description). In the selected observations, 
the f.o.v. always contains three sources: A0535, 
the Crab Nebula, and the low-mass X-ray binary 
4U 0614-1-091. However, the latter is much fainter 
than the former two, so we can safely neglect it 
from our sky model. For each revolution, we ex¬ 
tract the background-subtracted spectra of A0535 
and of the Crab, together with the correspond¬ 
ing response matrices with spidaE We selected 
spectral counts in the 20-135 keV range, and di¬ 
vided them into 45 logarithmically-spaced bins. 
For the evaluation of the spectral counts from each 
source we assume a variation time-scale of one sew 


^SPI Data Analysis Interface. The software has been devel- 
oped at the IRAP Toulouse, and it is available on the web 
at http:// sip:ina-2.ces r ,fr/i ntegral/spidai For more 
details, see lBurke et alj ll2014f) . 
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for the flux of A0535, while for the Crab we assume 
a constant flux over the time-scale of the observa¬ 
tion. We exclude sew affected by high particle 
background rates, e.g. in case of solar flares or 
close to the satellite ingress/egress from the Earth 
radiation belts. 

Thanks to their long exposures, we are able to 
split rev. 1019 in three sub-sets (labeled A, B, and 
C, respectively), and rev. 1021 in two sub-sets (la¬ 
beled respectively A and B), in order to check for 
variations of the spectral and timing parameters 
within the revolution time-scale. 

We also extract event lists corresponding to 
each data sub-set in order to follow the evolu¬ 
tion of the pulsar spin period and of the associated 
pulse profiles. 


3. RESULTS 


3.1. Spectral Analysis 


The analysis of the source spectra has bee n per¬ 
formed with XSPEC 12.8.2 ( Arnaudlll99fill . The 
hard X-ray emission of accreting pulsars is thought 
to originate from the thermal and bulk comp- 
tonization of seed photons in the a ccretion col¬ 
umn (see e.g. iBecker fc Wolfil I2007L for details). 
The resulting spectral shape can be well approx¬ 
imated by a power law continuum with exponen¬ 
tial cut-off at high energies (cutof fpl in XSPEC), 
F{E) = AE-'^ eyiY>{-E/Ecut), where P, E^ut, 
and A are respectively the spectral index of the 
power law, the folding energy and the normaliza¬ 
tion constant of the model spectrum, and E is the 
energy of the incoming photons in keV. We first 
perform a fit of each spectrum with this model. 
We notice that the spectrum of rev. 1028 has low 
count statistics, hence we exclude it from the spec¬ 
tral analysis. 


All spectra show large structured residuals 
around ^ 45keV (see Figured]), pointing out the 
presence of the fundamental cyclotron feature in 
our data. We account for this feature by includ¬ 
ing a multiplicative absorption line with Gaussian 
profile (gabs in XSPEC) in the spectral model, 
which becomes E'{E) = E'(E) exp[— t(E)], where 
t{E) = {TcyclV^cFcyc)e^I>[-{E - Ecy})"^/2aly^] 
and Ecyc, cTcyc, and Tcyc are respectively the cen¬ 
troid energy, the width, and the depth of the 
cyclotron line. This results in significant improve¬ 
ments in the ia all cases. 


We notice that there is a strong coupling be¬ 
tween the spectral index and the cut-off energy 
of the continuum model (Figured]). In addition, 
the width of the absorption line is anti-correlated 
with parameters of the continuum. This may re¬ 
sult in unrealistic values of the fit parameters. For 
this reason, we fix the s pectral index to an aver- 
age value, T = 0.5 (e.g. Miiller et al. 2013a^ and 
repeat the fits. 

To assess the significance of the absorption lines 
we proceeded as follows. First, for each data-set 
we generate 10000 simulated spectra with the in¬ 
built fakeit command from XSPEC, using the 
best-fit model as a template. Then we fit each 
simulated spectrum with either the cutoffpl or 
the cutoffpl * gabs models, imposing that the 
energy of the line is equal to the correspond¬ 
ing best-fit value, and taking note of the result¬ 
ing x^- From these we build the distribution of 
the F statistics, defined a s F = (Xo/^o)/(Xi/^i) 
(see Orlandini et al. 20121 and references therein), 
where i>i is the number of degrees of freedom of 
model i and the subscripts i = 0,1 correspond 
respectively, to the models without and with the 
absorption line. Finally we infer the probability 
of a chance improvement of the x^ by counting 
how many times the simulated values of F were 
larger than that obtained from the real data. In 
all cases the estimated chance probability is lower 
than 10“^, implying a significance greater than 
3.9cr. 

From part B of rev. 1021 onward the spectrum 
exhibits structured residuals at ^ 100 keV. The 
existence of a second cyclotron line at ^ 100 keV 
in the spectrum of this source ha s been reported 


sever al times in the past (see e.g. ICaballero et al 


2 OO 7 II . and the residuals observed in our spectra 
are likely due to this feature. As in the case of 
the line at lower energy, we perform simulations 
in order to assess the significance of these absorp¬ 
tion lines, where we now include the fundamental 
feature at ^ 45keV in the null hypothesis. The 
line is formally significant (> 3cr) for rev. 1021B, 
1022, 1023, and 1026, while for rev. 1024 the sig¬ 
nificance is barely ^ 2.3cr. We nevertheless include 
the second line in the spectral model for all these 
observations. Finally, for rev. 1025 and 1027 we 
find no evidence of residuals around 100 keV, pos¬ 
sibly due to a lack of statistic caused by the com¬ 
bined effect of the fading flux of the source and of 
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the short exposure times. We report the results of 
the spectral fitting in Table [21 where all the un¬ 
certainties are at the la confidence level and have 
been estimated by means of simulations. 


For each spectrum, we estimate the corre¬ 
sponding flux in the 20 — 80 keV range, and 
normalize it to that of the Crab Nebula ob¬ 
tained simultaneously from the same observa¬ 
tions (Figure |T]). The flux reaches a maximum 
of - 3.87X 10-®ergcm-2s-i (- 2.6 Crab) during 
rev. 1021B. This corresponds, at a distance of 
— 2 kpc, to a luminosity Lboi — 4.9 x lO^^ergs"^, 
where we assumed isotropic emissi on and a bolo- 
metric correction factor of 40% (Bildsten et al 


199711 . 


We then look for possible correlations of the 
spectral parameters with the observed flux, in or¬ 
der to constrain the accretion regime of the source. 
We find that the cut-off energy of the spectral con¬ 
tinuum is anti-correlated with the flux (Figured 
upper panel). Accordingly, the spectrum softens 
as the flux increases, as showed by the hardness 
ratio between fluxes in the 20 — 25 and 70 — 80 
keV bands (lower panel of Figure [4]). The bound¬ 
aries of the two energy bands were chosen in order 
to limit as much as possible the contribution of 
the two cyclotron lines. 


Intriguingly, the relations between the param¬ 
eters of the fundamental cyclotron line and the 
flux of the source differ between each other (Fig¬ 
ure HI left panels). The centroid energy Acyci is 
positively-correlated with the flux of A0535, going 
from — 50 keV at high luminosity (rev. 1021B) 
to — 44 keV at low luminosity (rev. 1026). We 
notice that the energy of the line at the lowest lu¬ 
minosity (rev. 1027) is compatible with that from 
the previous revolution, which implies a possible 
change in the correlation at about this luminosity. 


On the other hand, the width of the line (Tcyci 
remains in the 9—13 keV range, and it does not 
exhibit a correlation with the flux, while the line 
depth Tcyci presents a clear anti-correlation with 
the flux (see Figure [5]), going from — 13 at low 
luminosity to — 3 at high luminosity. Figure |6| 
shows the contour plot for the optical depth Tcya 
versus the centroid energy Ecya of the fundamen¬ 
tal cyclotron line at the beginning (rev. 1019A), 
maximum (rev. 1021A and 1022), and towards 
the end of the outburst (rev. 1026). It can be 
observed that the depth of the line exhibits sig¬ 


nificant variations during the ontburst, while the 
variations of the centroid energy are less signifi¬ 
cant. The parameters of the line at the end of the 
outburst are very close to the ones at the begin¬ 
ning with similar flux level, suggesting the absence 
of hysteresis effects during the outburst. 

Finally, we find no correlation between the flux 
of the source and the parameters of the first har¬ 
monic. However, the energy of the line seems to 
increase from — 100 keV to — 108 keV between 
rev. 1021B and 1022, and to decrease again to 
— 100 keV from rev. 1024 onward even though 
the lack of statistics prevents any firm conclusion. 


3.2. Timing analysis 


For each event list, we first convert photon ar¬ 
rival times to the Solar System barycenter and ap¬ 
ply corrections to account for the binary motion of 
A053 5. Ephemerid s for t he binary orbit are taken 
from [Finger et al.l (1996), adopting the updated 
values from Camero-Arranz et al.l ( 2012h for the 
epoch of the periastron passage and the orbital 
period. 

W e then carry out a 




test ( Buccheri et al 


19831) around the known pulse frequency of the 
source on each event list. For this task, we se¬ 
lect only photons in the 24 — 49 keV range, where 
the signal-to-noise ratio is higher. To further in¬ 
crease the signal-to-noise ratio, we use the infor¬ 
mation from the spectral deconvolution algorithm: 
for all the available sew, we estimate the fraction 
a of source photons that reach each detector due 
to the shadow casted by the coded mask on the 
detectors, and select only photons f rom detectors 


with a larger than a given value Isee lMolkov et al 


20101 for details). In this work, we exclude events 
from detectors with a < 0.2. 

The periodogram exhibits a very large peak 
around the reported spin period of the pulsar in all 
data-sets, with values of the power ranging from 
— 6 X 10^ to — 1.5 X 1 0"^ against an expected av¬ 
erage of 4. Following iLeahvI (1983), we fit the 
periodograms around the peak of the distri¬ 
bution with a Gaussian model in order to better 
estimate the spin frequency. The la uncertainties 
on the best-fit freque ncies were infe r red by adapt¬ 
ing equation B17 of [Ransom et al.l (|2002l i to the 
Z 2 case. 

Then, we perform phase-connected timing in 
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order to determine the spin evolution of the pulsar 
during the outburst. The pulse phase at a time t 
can be expressed as: 


(j){t) — (f>o + fo{t — to) + 2fo{t ~ (1) 


where / is the spin frequency, the dots represent 
derivation with respect to the time, and the sub¬ 
script 0 denotes quantities evaluated at the refer¬ 
ence epoch. We take MJD = 55616.202, the mid¬ 
point of rev. 1021A, as the reference epoch tg. 
Starting from a simple constant frequency model, 
f{t) = /o, for each data-set, we produce pulse pro- 
hles using events in the 27 — 49keV range, and 
measure the phase shift as a function of time with 
respect to the reference epoch. The phase is mea¬ 
sured at the minimum of the pulse profile. We find 
a reasonable fit using a polynomial model which 
includes the derivatives of the frequency up to the 
fourth order. The best-fit coefhcents of the timing 
solution are reported in Table [S] 


We can now infer the evolution of the spin fre¬ 
quency and its first derivative along the outburst 
(see respectively upper and lower panels of Figure 
El). The increase in the spin frequency induced 
during the accretion episode is evident. More¬ 
over, the pulsar was already spinning up during 
the rising part of the outburst (rev. 1019A to 
C), and it reached a maximum spin-up rate of 
(6.48±0.05) X 10-12Hzs-i at MJD=55617.25, i.e. 
about five days after the passage at the periastron. 


The rate of spin-up correlates with the flux 
of the source (Figure IH). This behavior is ex¬ 
pected in the case of accretion from a prograde 
accretion disk le.g. lGhosh fc Laniblll979ll . because 
both the spin frequency derivative rate and the lu¬ 
minosity are proportional to the accretion rate. 
Fitting the / — Flux( 20 - 80 keV) relation with a 
power law model, / ^ (Flux/Flux^ax)^) we ob¬ 
tain 7 = 1.24 ± 0.05, higher than the standard 
value of 6/7. 

Finally, we build pulse profiles for each data-set 
in four different energy bands, 20 — 27, 27 — 36, 
36 — 49, and 70 — 85 keV, respectively, by folding 
the data according to our estimated timing solu¬ 
tion. The profiles are then normalized to their av¬ 
erage value. We avoid the energy interval between 
49 and 70 keV because of t he presence of a strong 


background feature (see IWeidenspointner et al 


20031) which significantly reduces the signal-to- 
noise ratio of the pulsation. The resulting pulse 
profiles are shown in Figures|9]toE2l At high lumi¬ 
nosity, the pulse profile is made of two sub-pulses. 
The width of each sub-pulse is Acf) ^ 0.3, where 4> 
is the rotational phase in units of the spin period. 
The relative strength of the sub-pulses is roughly 
equal in the lowest energy range, 20 — 27 keV. At 
higher energies the second peak becomes fainter 
than the first one, which dominates the pulse pro¬ 
file above 49 keV. As the luminosity decreases, 
both sub-pulses weaken. However, it can be no¬ 
ticed that the change in the first sub-pulse is much 
more dramatic than in the second one. At ener¬ 
gies below 49 keV the profile shows a broad, single 
pulse (A/) ^ 0.8), while in the 70 — 85 keV range 
the lack of statistics does not allow to determine 
the presence of pulsations. 


4. DISCUSSION 

We analyzed several INTEGRAL/SPl observa¬ 
tions of the transient Be/X-ray binary A0535, per¬ 
formed during a type H outburst of the source oc¬ 
curred in February 2011. 

Based on the data taken during a previous 
outburst. iMiiller et al. ( 2013al) suggested that the 
lack of significant flux-related variations of the cy¬ 
clotron line energy in A0535 result from a peculiar 
accretion geometry, or from a peculiar alignment 
between the pulsar and an observer at the Earth. 
As an alternative, they suggested that A0535 may 
be in a different accretion regime with respect to 
other sources, where the luminosity is close or 
below the so-called Coulomb stopping limit (see 


Becker et al.l 120121 for details), and in any case 
below the critical value Lcrit ^ 6.8 x 10^^ erg s”^, 
estimated for this source at all except the bright¬ 
est outbu rsts of the sourc e regis tered so far. The 
results of ICaballero et al.l (j2012h for the outburst 
of February 2011 seem to support this picture. 
However, the results of our spectral analysis are in 
striking contrast with those reported above. We 
found a significant increase in the energy on the 
fundamental cyclotron line, which correlates with 
the flux of the source. This discrepancy might be 
due to the different energy range used in the spec¬ 
tral analysis, and/or by the different hypotheses in 
the modelization of the continuum emission, which 
might have a considerable effect on the measured 
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energy of the cyclot r on line s, see e.g. the discus¬ 


sion of Miiller et al. ( 2013bh . 


In any case, the obseryation of a positiye cor¬ 
relation between the energy of the fundamental 
cyclotron line and the flux confirms that A0535 
is accreting at sub-critical regimes. Interestingly, 
the correlation seems to fall off at the end of the 
outburst. This behayiour is expected when the ac¬ 
cretion flow passes from a Coulomb-stopping de- 
celeration reg i me a gas shock deceleration regime 


(IBecker et al.l 12012!) . In 


tion occurs between rey. 
is between ^ 1.2 x 10^^ 


our case, the transi- 
1026 and 1027, that 
ergs“^ and ^ 6.8 x 


10 ^®ergs in broad agreement with the yalue of 
7.4 X 10^^ e rgs~^ obtained from Equation 54 of 


Becker et al.l (|2012!) , where we assumed for A0535 
a magnetic field of 4 x 10^^ G. 

The spectral analysis also highlighted a soft¬ 
ening of the continuum emission as the flux of 
the source increased. A similar behayior of the 
hard (20 — 100 keV) X-ray spec trum was also re¬ 
ported bv iBildsten et al. (1997) from BATSE ob- 
seryations of the giant outburst of 1994. This 
is apparently in contrast with the hardening of 
the spe ctrum with flux reported by iMiiller et al. 
( 2013al) . However, the hardness ratio estimated 
by these authors were obtained from different en¬ 
ergy bands, which included the soft (5-20 keV) 
X-ray contribution as well. A global hardening of 
the spectrum could result from a shift of the emis¬ 
sion peak towards higher energies. This would in 
turn cause the softening of the spectrum observed 
above 20 keV. 

The spin history of A0535 has been tracked 
since its discovery fsee lCamero-Arranz et al.]l2012L 
and references therein). The pulsar shows a steady 
spin-down during quiescence, which is overlapped 
with spin-up episodes during giant outburst. From 
the timing analysis we obtained a phase-connected 
timing solution for the evolution of the spin fre¬ 
quency relative to the February 2011 outburst. 
Making use of this timing solution, we inferred 
the evolution of the spin-up rate, which reached a 
maximum value of / = (6.48±0.05) x 10“^^ Hz s“^ 
at MJD 55617.25, in good agreement with the 
value of / = (6.35 ± 0 .0 5) x 1 0~^^ Hz s~^ found by 


Camero-Arranz et al.l ( 2012h using Eermf/CBM 
data of the same outburst. We showed that the 
spin-up rate of A0535 during the outburst was cor¬ 
related with the flux of the source, / ^ (Flux)^'^^, 


and it is clearly different from the theoretical value 
of 6/7 expected in the case of disk accretion. 
Th is discrepancy h as be en recently reported also 
bv ISugizaki et al. ( 2015 ) using Eermi/GBM and 
MAXl/GSC data of A0535 and other Be X-ray bi¬ 
naries, and it could be explained by unaccounted 
variations of the beaming geometry and/or of the 
bolometric correction factor throughout the out¬ 
burst, which may affect the determination of the 
true luminosity of the source. 

Our study of the energy-dependent pulse pro¬ 
files and their evolution along the outburst gave 
results in line with thos e found in previous 


outbursts of the source (Bildsten et al.l 11997 


Camero-Arranz et al.l l2012l) . Th e pulse profiles 


of A0 535 have been analyzed by ICaballero et al 


(I 2 OIII) . us ing the beam patt ern decomposition 
method of Kraus et'^] ( 19951) . and their results 
are consistent with the emission from a hollow 
cone (i.e. the accretion column) above the polar 
cap(s) of the neutron star. The transition from 
a top-hat, energy independent profile at low lu¬ 
minosity to a highly asymmetric double-peaked 
and energy-dependent profile at high luminosity 
might be explained in terms, for instance, of vari¬ 
ation of the height and shape of the accretion 
column. The broad, top-hat and almost energy- 
independent profile at low luminosity suggests 
that the emission region occupies a large part 
of the neutron star surface (even accounting for 
light-bending effects) and that the physical condi¬ 
tions, and thus the optical depth, do not change 
much with the rotational phase. On the other 
hand, the asymmetric and energy-dependent pro¬ 
file observed at high luminosity might result from 
a non-symmetric accretion column, due to either 
offset dipole or non-dipolar magnetic fields, and 
thus, to a non-trivial and energy-dependent rela¬ 
tion between the physical conditions and geometry 
of the emission region and the rotational phase. 
This relation can be verified by means of phase- 
resolved spectroscopy, which will be detailed in a 
forthcoming paper (Sartore et ah, in preparation). 
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IJD (= MJD - 51544) 


Fig. 1.— Evolution of the of 20-80 keV flux of 
A0535 during the outburst of February 2011. 



100 100 100 
Energy [keV] Energy [keV] Energy [keV] 


Fig. 2.— Spectra (top panel) and associated fit 
residuals (second to fourth panel from top to bot¬ 
tom) of A0535 at different epochs of the out¬ 
burst, and corresponding, from left to right, to 
rev. 1019A, rev.1022, and rev. 1026, respectively. 
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Spectral Index 

Fig. 3.— Contour plot of the cut-off energy Ecut 
versus the spectral index F. 



Flux [Crob] 

Fig. 4.— ( Top) Relation between the best-fit cut¬ 
off energy of the continuum and the 20-80 keV 
flux. (Botton) Relation between the hardness ratio 
of the spectrum and the flux of the source. The 
soft and hard fluxes have been estimated in the 
20-25, and 70-80 keV bands, respectively. Error 
bars represent the Icr uncertainties. 
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Fig. 5.— Relation between the best-fit parameters 
and the 20-80 keV flux of A0535 for the cyclotron 
lines at ~ 45 keV (left), and ~ 100 keV (right). 
The top, middle and bottom panels show respec¬ 
tively the centroid energy, the width and the depth 
of the lines. Error bars represent the Icr uncertain¬ 
ties. 



Fig. 6 .— Contour plot of the line depth Tcyc^ 
versus the centroid energy Ecyc^ of the funda¬ 
mental cyclotron line. Countour levels are drawn 
at xLn + 2.30 (solid), xl,in + 4-61 (dashed), 
Xmin + 9.21 (dotted), respectively, and represent 
the Icr, 2(7, and 3 ct confidence intervals for the two 
parameters of interest. 



MJC [xl0‘ doys] 


Fig. 7.— Evolution of the spin frequency {top) 
and of the spin frequency derivative {bottom) of 
A0535 during the outburst. Data points in the up¬ 
per panel show the periods measured with the 
test. The solid lines represents the best-fit model 
obtained from timing analysis. The shaded area 
denotes the Ict error region. The vertical dotted 
line shows the epoch of the passage at the pe- 
riastron. The dashed vertical line in the lower 
panel shows the epoch corresponding to the max¬ 
imum spin frequency derivative (i.e. of maximum 
torque). 



Fluxpp _ 50 kev) [xl 0 '‘ erg s'’ cm'’] 


Fig. 8 .— Relation between the spin frequency 
derivative and the 20-80 keV flux of A0535. 
Overplotted are two power law models, / 
(Flux/Fluxmax)"'', with 7 = 6/7 (dotted) and 7 = 
1.24 (dashed), respectivley. 
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Fig. 9.— Energy-dependent pulse profiles of A0535, normalized to their average value. 
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Fig. 10.— ...continued from Figure[9] 
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Fig. 11.— ...continued from Figure [TUI 
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Fig. 12.— ...continued from Figure fTT] 
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Table 1 

Log of integral observation of A0335+26 during the outburst of February 2011 


Revolution 

Start date of the observation 
[YYYY - MM - DD] 

Total exposure 
[ks] 

1019 

2011-02-16 

186.16 

1021 

2011-02-23 

101.71 

1022 

2011-02-27 

53.50 

1023 

2011-03-02 

43.12 

1024 

2011-03-05 

43.88 

1025 

2011-03-08 

46.14 

1026 

2011-03-11 

64.33 

1027 

2011-03-14 

50.33 

1028 

2011-03-17 

50.49 
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Table 2 

Best-fit spectral parameters for the INTEGRAL/SVl observations of A0535. 


Data-set 

1019A 

1019B 

1019C 

1021A 

1021B 

1022 

1023 

1024 

1025 

1026 

1027 

r 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

Ecut 


16.4t«;i 

l6.8«;i 

l5.7l«;i 

l5.8l«;= 

l5.6l«;l 

l5.9l«;l 

l6.5l«;= 

16.4l«;= 

i7.5iS:2 

l6.8«;t 

Ecyci 

45.6t°;J 

45.6lS:7 

47.4i°;® 

48.2;«;| 

49.8lJ;« 

46.4l«;« 

46.81“;^ 

47.9to.l 

46.lj;i'° 

43.9lS:t 

44.9to.^ 

^cyci 


8.9«;® 


9.8tS:I 

11 1+1.5 

o O+0.7 

o-o_0.6 

9.41°;? 



y-o_0.6 


Tcyci 

IS.ltf.o 

4 9+0-8 

« 0 + 1-2 
0-0_Q g 

o ^-(-0.5 
'^•0-0.4 

4 

^•"^-0.6 

9 8+0-4 

3 9+0-0 


8-011.3 

lo-alj'i 

ll-dti.g 

Significance 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

> 3.9(7 

Ecyc2 

- 

- 

- 

- 

98-611.I 

107.11^;° 

105 . 4 ;^ = 

100.3tt,® 

- 

99.2tJ;| 

- 

^cyc2 

- 

- 

- 

- 


P 0 + 1-5 

0.0_i 7 


<7 0 + 3.8 
' •^-3.2 

- 



Tcyc2 

- 

- 

- 

- 

l7.0i^;^ 

23.0j:"/ 

29.ol:^Y 


- 

r\A 7+896.1 


Significance 

- 

- 

- 

- 

3.10- 

3.6(7 

> 3.9(7 

2.3(7 

- 

3.5(7 

- 

Flux 

[Crab] 

1.14 

1.40 

1.76 

2.49 

2.62 

2.46 

2.20 

1.71 

1.21 

0.80 

0.36 

Xr\d(d.o.f.) 

2.34(37) 

1.84(37) 

2.84(37) 

1.20(36) 

1.09(25) 

1.08(34) 

1.41(34) 

1.12(34) 

1.80(37) 

1.50(34) 

1.63(27) 



Table 3 

Best-fit parameters of the polynomial fit describing the evolution of the spin 
PARAMETERS OF A0535 DURING THE OUTBURST. QUOTED ERRORS ARE THE FORMAL la 

UNCERTAINTIES FROM THE BEST-FIT. 


Epoch 

[MJD] 

55616.202 

fo 

(9.6793 ± 0.0001) X 10“^ Hz 

fo 

(6.43 ± 0.05) X 10“^^ Hzs“^ 

fo 

(1.21 ± 0.17) X 10“^® Hzs“^ 

fo 

(-1.43 ± 0.09) X 10“^® Hzs-® 

Jo 

(1.67±0.15) X 10"®® Hzs-* 
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